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1 | INTRODUCTION 


Summary 

Porcine deltacoronavirus (PDCoV) is a newly discovered coronavirus that causes 
diarrhoea, vomiting and dehydration in sucking and nursing piglets. It was first 
reported in Hong Kong in 2012 and has since been discovered in the United States, 
Canada, South Korea, mainland China, Thailand and Laos. PDCoV has been experi- 
mentally proved to lead to diarrhoea in swine and it was detected positive in pigs in 
Guangdong, southern China. In our study, 252 faecal and intestinal samples from 
sucking piglets and sows with diarrhoea were surveyed for common enteric viruses. 
We found a prevalence of PDCoV (21.8%), porcine epidemic diarrhoea virus 
(65.5%), transmissible gastroenteritis virus (0%), rotavirus group A (25.0%) and por- 
cine kobuvirus (68.7%). We isolated 13 PDCoV strains and discovered that PDCoV 
infections were often co-infections with kobuvirus rather than the commonly linked 
porcine epidemic diarrhoea virus. Phylogenetic analysis of S gene and N gene 
revealed that 11 of 13 PDCoV strains belonged to Chinese lineage. As for the left 
two strains, one single strain (CHN-GD16-05) belonged to American and Korean lin- 
eages while another strain (CHN-GD16-03) was similar to a Thai strain, but only in 
the S gene. This suggested a possible recombination event between the Thai and 


the newly described Chinese strain. 
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(PHEV) is a Betacoronavirus and porcine deltacoronavirus (PDCoV), a 


Deltacoronavirus. 


Coronaviruses belong to Coronavirinae subfamily within Coronaviri- 
dae family and the Nidovirales order. They are spherical, enveloped 
RNA viruses that cause enteric, respiratory and other diseases in 
humans and animals. In 2012 after the discovery of the Deltacoron- 
avirus in animal species in Hong Kong (Woo et al., 2012), the classi- 
fication of Coronavirinae subfamily was expanded. There are 
presently four genera in this subfamily: Alphacoronavirus, Betacoron- 
avirus, Gamacoronavirus and _ Deltacoronavirus. Five coronavirus 
strains have been isolated from pigs. These include the Alphacoron- 
avirus species porcine epidemic diarrhoea virus (PEDV), transmissible 
gastroenteritis virus (TGEV) and porcine respiratory coronavirus 


(PRCV). Porcine haemagglutinating encephalomyelitis coronavirus 


Porcine deltacoronavirus was first identified in swine specimens 
collected in 2009 in Hong Kong (Woo et al., 2012). It was identified 
later in several states of USA (Li et al., 2014; Marthaler, Jiang, Col- 
lins, & Rossow, 2014; Wang, Byrum, & Zhang, 2014) and has now 
spread to South Korea, mainland China, Thailand and Laos (Lee & 
Lee, 2014; Lee et al., 2016; Lorsirigool et al., 2016; Madapong et al., 
2016; Song et al., 2015). The emergence of PDCoV was most often 
associated with clinical infection symptoms including diarrhoea, vom- 
iting and dehydration in pigs, especially in sucking and nursing pig- 
lets (Janetanakit et al., 2016; Ma et al., 2015; Marthaler, Raymond 
et al. 2014; Song et al., 2015). These worldwide outbreaks of 


PDCoV have caused economic losses to the global commercial swine 
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industry (Chen, Gauger et al., 2015; Chen, Zhu et al., 2015; Hom- 
wong et al., 2016; Song et al., 2015; Wang et al., 2014). 

RT-PCR detection methods have provided data for the preva- 
lence of PDCoV in China (Dong et al., 2015; Song et al., 2015; Zhai 
et al., 2016). A previous study (Zhai et al., 2016) reported that the 
current prevalence in Southern China was 1.54% (5/390), and three 
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of the five detected strains were from Guangdong Province, reveal- 
ing that the information regarding the molecular epidemiology of 
PDCoVs in Guangdong was still limited. The current study further 
investigates the prevalence, epidemiology and genomic properties of 
PDCoV in Guangdong. 


2 | MATERIALS AND METHODS 


2.1 | Sample collection and detection methods 


In order to monitor the prevalence and sequence properties of 
PDCoV in Guangdong Province, 252 faecal and intestinal samples 
were collected from 11 commercial swine farms between December 
2015 and June 2016. These samples were all collected from sucking 
piglets and sows. We screened for PDCoV, PEDV, TGEV, Rotavirus 
A and porcine kobuvirus (PKV). 

Thirteen PDCoV cell culture isolates were propagated in swine 
testis (ST) cells and plaque purified twice. These isolates were 
obtained from case submissions from the 252 swine samples we col- 
lected previously. Total RNA was extracted from samples using a 
commercial kit following the manufacturer’s recommendations (Axy- 
gen Scientific). All RNA samples were stored at —80°C. RNA was 
reverse-transcribed using a PrimeScript RT Reagent Kit with gDNA 
Eraser (Takara, Biotechnology, Dalian, China) following the manufac- 
turer's instructions. PCR primers (5’-ATCCTCCAAGGAGGCTATGC-3’ 
and 5’-GCGAATTCTGGATCGTTIGTT-3’) were designed for the viral 
M gene of PDCoV strain CHJXNI2/2015/China (GenBank accession 
KR131621) (Song etal., 2015). The thermal cycling (worked on 
Bio-Rad T100, Forster City, CA) was operated with the following 
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thermal profile: 94°C for 5 min, 30 cycles of 94°C for 45 s, 55°C for 
30 s, 72°C for 1 min and a final step of 72°C for 10 min. Amplicons 
were analysed on 1% agarose gels. Moreover, on the basis of the 
four porcine enteric pathogens below, additional RT-PCR detections 
from the 13 tissue cultured purified PDCoV-positive samples target- 
ing porcine bocavirus (PBV), porcine sapelovirus (PSV) and porcine 
astrovirus (PAstV) were added. Primers specific for PEDV, TGEV, 
Rotavirus A, PKV, PBV, PSV and PAstV have been described at 
Table S1. 


2.2 | Amplification, cloning and sequencing the 
Spike (S) protein and Nucleocapsid protein (N) genes 


The complete S and N genes were amplified using three pairs of pri- 
mers: PDCoV-S1, 5’- ATGCAGAGAGCTCTATTG-3’ and 5’-TATTT 
CAACTTCGCCATC-3’; PDCoV-S2, 5’-CGACCATCCATAGTTTCA-3’ 
and PDCoV-S2, 5’-CTACCATTCCTTAAACTT-3’; and PDCoV-N, 5’- 
ATGGCCGCACCAGTAGTC-3’ and 5’-CTACGCTGCTGATTCCTG-3’ 
(based on CHJXNI2/2015/China). PCR amplification was carried out 
using the LA Taq polymerase kit (Takara, Biotechnology, Dalian, 
China) using directions supplied by the manufacturer. PCR products 
were purified using the Gel Band Purification Kit (Omega Bio-Tek, 
USA) and then cloned into the PMD-19T vector (Takara) using an 
In-fusion PCR Cloning Kit (Takara). The recombinant plasmids were 
Institute (Shenzhen, 


sequenced by the Beijing Genomics 


Guangdong). 


2.3. | Sequence alignment and phylogenetic analysis 


Nucleotide sequences were submitted to GenBank with accession 
numbers (KYO78891—KY078903 for nucleocapsid gene and 
KY078904—-KY078916 for spike gene of CHN-GD16-01 to CHN- 
GD16-13, respectively). Reference sequences included 33 strains for 
the S gene and 33 strains of N gene from different farms of global 
isolates obtained from GenBank (Table S2 and S3). DNA sequences 


TABLE 1 RT-PCR detection of porcine enteric viruses PEDV, PDCoV, Rotavirus group A, PKV and TGEV in diarrhoeal faecal or intestinal 


samples from December 2015 to June 2016 in Guangdong Province 


Diarrhoeal samples no. (%) PDCoV no. (%) 


Months of year SP/S Total SP S 


December, 2015 36/17 (67.9/32.1) 53 (21.0) 10 4 
January, 2016 26/18 (591/401) 4417.5) 9 6 
February, 2016 31/5 (86.1/13.9) 36(14.3) 7 @ 
March, 2016 dy 1(9S:e/6.2) Weyl (exGisy 2 1 
April, 2016 34/2 (94.4/55.6)  36(14.3) 2 0 
May, 2016 t3/3 78.0) 2107) zo (9) 5 0 
June, 2016 43/1 (97.7/2.3) 44(17.5) 9 0 
Total 203/49 (80.6/19.4) 252 44 (21.7) 11 (22.4) 
BS) (alse) 


Rotavirus A no. 


PEDV no. (%) (%) PKV no. (%) 

SP s SP S SP S TGEV 

22 9 12 2 27 6 @ 
14 ali 4 @ is ddl @ 
24 3 0 28 5 0 
‘Lal 0 4 @ 13 2 @ 
27 il 13 0 20 il 0 
a il 4 © cal S © 
32 0 16 0 26 5 0 
39 (68:5) 2653.1) rol (300)) 2 (4.1) al40(6910)= 3367.3)" 0 

165 (65.5) 63 (25.0) 7 3468.7) 


PEDV, porcine epidemic diarrhoea virus; PDCoV, porcine deltacoronavirus; PKV, porcine kobuvirus; TGEV, transmissible gastroenteritis virus; SP, 


Sucking piglets (<20 days); S, Sow. 
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were assembled using the Clustal W program of DNAStar V7.1 
(Madison, WI, USA). Phylogenetic trees were constructed using the 
neighbour-joining method in MEGA 6.0 software with bootstrap 
analysis of 1,000 replicates. The percentage of replicate trees in 
which the associated taxa clustered are shown as nearby branches 
(Tamura, Nei, & Kumar, 2004; Tamura, Stecher, Peterson, Filipski, & 
Kumar, 2013). 


2.4 | Prediction of epitopes for spike protein 
antigen sequence and N-glycosylation of PDCoV 


The secondary structure and surface properties of S1 protein of 
PDCoV were analysed with Protean, DNAStar V.7.1 (Madison, WI, 
USA) to predict its antigenic epitopes. The secondary structures of 
these S proteins were analysed by Garnier-Robson and Chou—Fasman 
methods, and surface properties, including hydrophilicity, flexible 
regions, accessibility and antigenicity were analysed by Kyte—Doolittle, 
Karplus—Schulz, Emini and Jameson—Wolf methods, respectively. Also, 
further prediction of epitope positions were confirmed on online ser- 
ver at http://www.cbs.dtu.dk/services/BepiPred/. Determination of 
potential N-glycosylation sites of the S protein was predicted with 


software available at http://www.cbs.dtu.dk/services/NetNGlyc. 


3 | RESULTS 


3.1 | PDCoV detection 


Samples (N = 252) from 11 swine farms in Guangdong Province 
were collected during a 7-month period. Approximately 22% were 
PDCoV positive while 65%, 25% and 69% contained PEDV, Rota A 
and PKV, respectively. Transmissible gastroenteritis virus was not 
detected in any samples. Sucking piglets showed slight lower PDCoV 
infection rate (21.7%, 44/203) than sows (22.4%, 11/49), while other 
types of pathogens are on the contrary (Table. 1). In addition, 
co-infections with PEDV and PKV were the most common (44.4%, 
79/178) (Table 2). PDCoV occurrence was examined further by an 
additional RT-PCR analysis using 13 PDCoV isolates we successfully 
purified in cell culture (Table 3). All 13 samples were collected from 
faeces and intestinal contents of sucking piglets. PDCoV-PEDV-PKV 
and PDCoV-PKV co-infection occurred with high frequency; Rota A 
was detected in three samples and PBV in only one. As for other 
porcine enteric pathogens, such as TGEV, PSV and PAstV, they were 


all examined negative in this work. 


3.2 | Homology and phylogenetic analysis of S gene 
and N gene of PDCoV isolates 


The 13 PDCoV genomes were designated CHN-GD16-01 to 
CHN-GD16-13. For S gene, sequencing results revealed that CHN- 
GD16-03 and CHN-GD16-05 possessed 3,483 nucleotides, 
CHN-GD16-06 was included 3,477 nucleotides and other 10 strains 
composed of 3,480 nucleotides in size. Interestingly, the viral strain 
CHN-GD16-06 contained unique AAT and GTT deletions compared 
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TABLE 2 Co-infections of different pathogens in samples 
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Virus PDCoV- 
positive %, positive %, 
samples no. (n = 246) no. (n = 55) 
One pathogen only 
PDCoV 14 20.6 14 
PEDV 32 47.1 
Rota A Z. D8: 
TGEV 0) 0) 
PKV 20 29.4 
Total 68 27.6 14 23.5 
2 pathogens 
PDCoV+PEDV 2 1.6 2 74 
PDCoV+Rota A 0) @) 0) 0) 
PDCoV+PKV 25 19.3 25 926 
PEDV+Rota A s) 23 
PEDV+PKV 7? 61.2 
Rota A+PKV 20 15:5 
Total 129 52.4 27 49.1 
3 pathogens 
PDCoV+PEDV+Rota A @) 0) 0) 0) 
PDCoV+PEDV+PKV 11 2.9 it 100 
PDCoV+Rota A+PKV 0) 0) 0) 0) 
PEDV+Rota A+PKV o5) Ton 
Total 46 13.7 11 83.6 
4 pathogens 
PDCoV+PEDV+Rota 3 1.2 3 a 


A+PKV 


PDCoV, porcine deltacoronavirus; PEDV, porcine epidemic diarrhoea 
virus; TGEV, transmissible gastroenteritis virus; PKV, porcine kobuvirus. 


to the US and individual Asian strains. Their nucleotide identities 
were represented in dispersed clusters within the reference strains 
(Table S2). From this point, we used only the last two numbers of 
the names of our cultivated strains to eliminate confusion with 
strains from other studies, that is, CHN-GD16-01 = strain 01. Most 
of the viral strains (01, 02, 09, 12 and 13) had a high nucleotide 
similarity (0.985-0.990) with CHJXNI2/China/2015 (Song et al., 
2015). On the other hand, 04, 06, 07, 08, 10, 11 shared 0.988— 
0.992 similarity with CHN/JS/2014 (Dong et al., 2015). Unlike the 
other strains from this work, O03 had the highest nucleotide homol- 
ogy (0.971) with HKU15-44, a strain first identified from swine in 
Hong Kong. Interestingly, strain O5 had a maximum _ nucleotide 
identity (0.985) with KNU14-04 (Lee & Lee, 2014) from South 
Korea. Moreover, the neighbour-joining phylogenetic tree recon- 
structed from the PDCoV S genes in isolate 03 was located at the 
Thai branch. The others were clustered within the Chinese branch. 
Strain O5 was located at a single branch next to the US and Korean 
branches (Figure 1). 

Within the N= gene, our 13 purified samples contained 
1,029 nt indicating a lack of insertions and deletions. These 13 
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TABLE 3 _ Detection of porcine enteric pathogens from 13 tissue cultured purified PDCoV-positive samples from 11 swine farms in 


Guangdong 

Farms Sample’s name Date of isolation Specimens _ Pig herds 

A GD16-01 24 February 2016 ~—Faeces Sucking Piglets 
B GD16-02 15 March 2016 Faeces Sucking Piglets 
C GD16-03 18 March 2016 Intestine Sucking Piglets 
D GD16-04 14 March 2016 Intestine Sucking Piglets 
E GD16-05 5 January 2016 Intestine Sucking Piglets 
FE GD16-06 18 February 2016 Intestine Sucking Piglets 
G GD16-07 23 January 2016 Intestine Sucking Piglets 
G GD16-08 30 January, 2016 Intestine Sucking Piglets 
H GD16-09 4 February 2016 Intestine Sucking Piglets 
| GD16-10 2 January, 2016 Intestine Sucking Piglets 
| GD16-11 4 February 2016 Intestine Sucking Piglets 
J GD16-12 2 January 2016 Intestine Sucking Piglets 
K GD16-13 26 January 2016 Intestine Sucking Piglets 


MAI ET AL. 
Porcine enteric pathogens 
PDCoV PEDV TGEV- RotaA PKV’ PBV_ PSV_ PAstV 
+ a a a a + — NT 
+ = = = + = = — 
+ + — + + = — — 
i + _ — +r = NT NT 
+ _ _ — + — NT NT 
+ + — _ +r = NT NT 
i. a = = + — NT NT 
+ + = + = — NT NT 
+ + — — 7 7 NT NT 
i +t _ _ + = NT NT 
ot + — — + 7 NT NT 
+ + - a +r = NT NT 
+ + — + + — NT NT 


PDCoV, porcine deltacoronavirus; PEDV, porcine epidemic diarrhoea virus; TGEV, transmissible gastroenteritis virus; PKV, porcine kobuvirus; PBV, por- 


cine bocavirus; PSV, porcine sapelovirus; PAstV, porcine astrovirus; +, positive; —, negative NT, not tested. 


Sucking piglets were all less than 20 days old. 


samples were also had the highest nucleotide identities with dif- 
ferent PDCoV strains. For example, strain 05 was most closely 
to KNU14-04 (S. Korea), OH11846 and Minnesota292 (0.993) 
(USA). The other strains were most closely to CH/SXD1/2015, 
Sichuan/S27 and CHN/GX01/2015 (Chen, Gauger et al., 2015; 
Chen, Zhu et al., 2015; Wang, Yue, Fang, & Huang, 2015; Zhai 
et al., 2016) originating in different regions in China (Table S3). 
The phylogenetic tree of the N gene indicated that strain O5 
was located at the US branch, while other positive samples from 
this work were all clustered in a group within the Chinese 
branch (Figure 2). 

Overall, combining the results of homology and phylogenetic 
analysis of the S and N genes, we concluded that most of the cur- 
rent PDCoV strains in Guangdong Province are closely related to 
other Chinese PDCoV strains. Only a few of PDCoV identities in 


strains O3 and 05 might have evolved from different ancestors. 


3.3. | Spike protein B-cell antigenic epitope and 
N-glycosylation site prediction 


We identified 11 potential and distinct B-cell antigenic epitopes in 
the S1 protein (Figure $1). Peptides less than 4 amino acid homology 
were not computed. These predicted epitopes all contained mostly B 
structure and high hydrophilicity, strong accessibility, good flexibility 
and a high, predicted antigenicity. We also identified 17 potential N- 
glycosylation sites in most of the isolates as well as the reference 
strains. Strain 02, USA/IL/2014 and Korea/SL5/2015 each contained 
16 glycosylation sites with the lack of an N-glycosylation site at 
position 162. An additional deletion at position 41 was observed in 
Thailand strains PDCoV/S5011/2015, PDCoV/S5015L/2015 and 
PDCoV/S5024/2015. 


4 | DISCUSSION 

Until now, numerous reports had revealed that PEDV was the most 
common porcine enteric pathogen in pigs due to its high prevalence 
in diarrhoeal samples. PDCoV-PEDV co-infection occurred at a rate 
greater than PDCoV-Rota A and PDCoV-TGEV. But in our study, 
co-infection of PDCoV and PKV was greatly higher than PDCoV 
and PEDV co-occurrence (45.4%, 25/55 vs. 3.6%, 2/55, Table. 2); 
additionally, PDCoV-PEDV-PKV co-infection also processed certain 
proportion on PDCoV-positive samples (20%, 11/55). PKV was also 
common pathogens causing diarrhoea and viral gastroenteritis in 
pigs, especially in piglets (Di, Angeloni, Tofani, Monini, & Ruggeri, 
2015; Pfankuche et al., 2016) and it showed a high co-infection 
rate with PEDV (Wang, Lan, & Yang, 2016; Zhao et al., 2016). PKV 
detection has not been included with most other studies of PDCoV; 
therefore, evidence is lacking to prove that PKV co-infects more 
frequently than PEDV with PDCoV. This warrants further 
investigation. 

HKU15-44 and HKU15-155 were the first two PDCoV strains 
that were sequenced (Woo etal., 2012). The whole S_ gene 
(3,483 nt) of published porcine deltacoronavirus isolates was anal- 
ysed by pairwise alignments. Most of the Chinese strains (HKU15- 
155, HB/2014, JS/2014, HN/2014, SXD1, CHJXIN2, GDO1/2015, 
GD02/2015, HN0O1/2015 and GX01/2015) contained a 3-nt (AAT) 
deletion (Chen, Gauger et al., 2015; Chen, Zhu et al., 2015; Dong 
et al., 2015; Song et al., 2015; Wang et al., 2015; Zhai et al., 2016). 
In strain CHN/AH/2004 (Dong et al., 2015), a 3-nt (TAA) deletion 
existed in its 3’-UTR. The N gene, however, lacked any nucleotide 
deletions or insertions. 

In our work, pairwise alignment analysis revealed that two S 


genes from our strains 03 and O5 lacked any deletions or insertions. 
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FIGURE 1 Phylogenetic analysis using the PDCoV spike protein (S) gene from different species. The S gene strains identified in this work 
are indicated with “black triangles”. GenBank accession numbers, countries and years of isolation are listed at the right side of each reference 
strain. The tree was constructed using MEGA 6.0 software with neighbour-joining (NJ) methods and 1,000 replicate sets on bootstrap analysis. 
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FIGURE 2 Phylogenetic analysis of the PDCoV N gene isolated from different species. The N genes from the strains of this work are 
indicated with “black diamonds”. GenBank accession numbers, countries and years of isolation are listed at the right side of each reference 
strains. The tree was constructed as per Figure 2, above. 
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The remaining 11 isolates all contained an AAT deletion resulting in 
an asparagine deletion. Furthermore, strain O06 possessed a GTT 
deletion, causing a valine deletion. More importantly, a C3447T sub- 
stitution of its nucleotide sequence changed a glutamine codon into 
a stop codon, resulting in an 11 amino acid (QPTPSFKFKEW) dele- 
tion. Although amino acid deletions and insertions have been 
reported, it remains unknown whether these have any biological sig- 
nificance to PDCoV pathogenesis. 

Two pairs of strains O7 and 08 from farm G and 10 and 11 from 
farm | were isolated from the same commercial swine farms, but at 
different times. Unlike O7 and 08 that were identical in nucleotide 
sequence, 10 and 11 showed nucleotide identities of 99.0% and 
98.3%, respectively, resulting in 10 amino acid changes. Interestingly, 
these differences were in the S protein and concentrated within the 
S1 region. According to the two phylogenetic trees established for 
the S and N genes, 11 of 13 strains clustered within the Chinese 
category and had the highest nucleotide identities with five pub- 
lished Chinese strains. This indicated that they were more closely 
related to other strains from China. 

Strain 05 isolated from farm E belonged to the US and Korean 
branch with the greatest identity with Korean strain KNU14-04 (Lee 
& Lee, 2014). As all reported PDCoV isolates from China are clus- 
tered into a single Chinese group, this is the first case of a US and 
Korean strain detected in diarrhoeal samples from pigs in China. 
Strain O3 is more related on a genomic scale to the Thai strains and 
had the lowest nucleotide identity with CHN/GDO2/2015 and 
CHN/GX01/2015 (95.8%, Table $2) (Zhai et al., 2016), which were 
collected from Guangdong and nearby Guangxi provinces. However, 
the N gene of this strain revealed a completely reverse conclusion; 
that it belongs to the Chinese group with a low nt identity (96.7%) 
with the Thai strains $5011 and 5015L (Janetanakit et al., 2016) 
(Fig. 2 and Table S3). Therefore, we inferred that a recombination 
event might occur in PDCoV strain 03. 

Recombination events are often reported in PEDV studies and the 
majority are intrarecombinants with different strain lineages of the 
same kinds of enteric coronaviruses (Jarvis et al., 2015; Tian et al., 
2014). A recent study (Boniotti et al., 2016; Valerij et al., 2016) identi- 
fied a swine enteric coronavirus (SeCoV) from Italy as more closely 
related to PEDV in the S gene, while the remainder of the genome 
shared highest identity with TGEV. Recombination events have not 
been reported in PDCoV strains. Since we were lack of the full-length 
genome of this PDCoV strain, the mechanism of the construction of 
this chimera as well as with the recombination in strain 03 is unknown, 
something suppose to be further work and explain with it. 

The spike protein is an important surface glycoprotein of coron- 
avirus and plays a significant role in receptor binding and fusion of 
the viral and cellular membranes (Chakraborti, Prabakaran, Xiao, & 
Dimitrov, 2005; Schwegmann-Wefels et al., 2009). And also, it 
mediates interspecies transmission (Bosch, Van, de Haan, & Rottier, 
2003). All coronavirus spike proteins share the same two functional 
components: an N-terminal subunit and a membrane-anchored sub- 
unit (C-terminal) that are covalently bound. The S protein of PDCoV 
can be similarly subdivided into $1 (1-573 aa) and S2 (574-1161 aa) 
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regions (Thachil, Gerber, Xiao, Huang, & Opriessnig, 2015). $1 is a 


dominant viral antigen and an ELISA is available for its detection 
(Thachil et al., 2015). Three different groups of S1 proteins from co- 
viral infections shared less than 12% amino acid sequence identities 
with each other (Wang, Deng et al., 2016). The B-cell response is 
directed against the spike protein of coronaviruses (Cao et al., 2015) 
and it plays an important role in pathogenesis of virus infection. 

In our study, 11 potential B-cell antigenic epitopes of PDCoV 
(Figure S1) and 16 or 17 N-glycosylation sites were analysed. Inter- 
estingly, in amino acids 549-561 in strain 05, a T5591 mutation 
greatly reduced the predicted level of its antigenicity, hydrophilicity, 
surface probability and flexibility (Figure S2 & Table S4). Whether 
this mutation has altered its antigenicity will be interesting to test in 
the laboratory. Overall, further studies are needed to confirm 
whether these alterations of B-cell antigenic epitopes and N-glyco- 
sylation sites affect the pathogenicity and antigenicity of each 
PDCoV strain. 

In conclusion, diarrhoeal samples collected from pigs in Guang- 
dong Province were screened to detect the prevalence of PDCoV. 
Phylogenetic analysis suggested that nearly all of the strains in main- 
land China were clustered into the Chinese lineage except one newly 
discovered PDCoV strain that had a close relationship with US and 
Korea strains. This complements the geographical lineage theory of 
global PDCoV distribution. The presence of another suspected 
recombinant strain will provide additional data to examine the diver- 


sity of the PDCoV genome. 
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